Quantum spin liquids and the metal-insulator transition in doped semiconductors 
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We describe a new possible route to the metal-insulator transition in doped semiconductors such 
as Si:P or Si:B. We explore the possibility that the loss of metallic transport occurs through Mott 
localization of electrons into a quantum spin liquid state with diffusive charge neutral "spinon" 
excitations. Such a quantum spin liquid state can appear as an intermediate phase between the 
metal and the Anderson-Mott insulator. An immediate testable consequence is the presence of 
metallic thermal conductivity at low temperature in the electrical insulator near the metal-insulator 
transition. Further we show that though the transition is second order the zero temperature residual 
electrical conductivity will jump as the transition is approached from the metallic side. However 
the electrical conductivity will have a non-monotonic temperature dependence that may complicate 
the extrapolation to zero temperature. Signatures in other experiments and some comparisons with 
existing data are made. 



Phenomena near the metal-insulator transition (MIT) 
in doped semiconductors such as Si:P or Si:B have been 
studied extensively for more than three decades [1-4]. 
Nevertheless, several aspects of the physics, for instance 
the detailed critical behavior [4-6], remain mysterious. In 
this paper we explore and develop the consequences of a 
new possible route to the MIT where a quantum spin 
liquid insulator appears as an intermediate phase be- 
tween the metal and the Anderson-Mott insulator. In 
recent years such a quantum spin liquid Mott insulator 
has been observed to intervene between the Fermi liq- 
uid metal and conventional magnetically ordered Mott 
insulators in a few different clean materials [7-9] . Here 
we study the strongly disordered situation appropriate 
to doped semiconductors and describe a variety of exper- 
imental consequences. 

When P is doped into Si, the extra electron of P forms 
a Hydrogen-like state with an effective Bohr radius of 
about a s=s 20 A[l]. A simple picture of the doped semi- 
conductor is as a collection of randomly placed "Hydro- 
gen" atoms. The system may then be described as a half- 
filled Hubbard model on a random lattice supplemented 
by the inclusion of the long range Coulomb interaction 
Vij between the electrons: 

H = - ^3 ( c L c ja + h.cj +U^2 n^n^ + ^VijUirij 

ij;a i ijtj 

(1) 

At low concentrations, the tij w toe Ti il a are small, 
the on-site U dominates and a Mott insulator of local 
moments results. The local-moments are coupled anti- 
ferromagnetically, and due to their random placement, 
preferentially form singlets with their closest available 
neighbor. The resulting random-singlet phase has an 
extremely broad distribution of singlet binding energies, 
giving rise to a diverging density of states for low-energy 
spin-excitations, contributing a divergent coefficient of 
heat capacity, 7 = C/T and spin-susceptibility x[10]- 
As the concentration of dopants, n, is increased, even- 



tually the typical Uj dominates over the U and a diffusive 
metal is obtained. A continuous phase transition between 
metal and insulator occurs at some critical intermediate 
concentration, n c , where Uj w U. Because of the random 
placement of dopants, a fraction of the local moments 
are very weakly coupled to the conducting electrons and 
survive unscreened into the metallic phase. The diffu- 
sive metal appears to be well described by a "two-fluid" 
model where the conducting electrons exist essentially de- 
coupled from a random fraction of weakly-coupled local- 
moments[l, 11]. As in the insulating phase, these lo- 
cal moments continue to dominate the low-temperature 
thermodynamic and magnetic properties of the metallic 
phase, but do not appear to strongly modify its transport 
properties. 

It is natural to ask: What is the fate of the con- 
ducting fluid across the metal-insulator transition? The 
conventional answer, implicitly adopted by most existing 
work[2, 3], is that all electron degrees of freedom are lo- 
calized by disorder[12], which is perturbativcly enhanced 
by interactions. In this scenario, shown in Fig. la, de- 
creasing n < n c gives a localized Anderson-Mott insu- 
lator with non-zero average density of states. As n is 
further decreased, the system crosses over continuously 
towards a correlation driven Mott-insulator of local mo- 
ments. 

In this paper, we point out a new and conceptually 
distinct scenario for the metal-insulator transition in 
doped semiconductors. In this scenario, the charged- 
conducting fluid is localized into a gapless quantum spin- 
liquid, but the electron thermal transport remains diffu- 
sive into the weakly insulating regime. There is grow- 
ing theoretical and experimental evidence that such gap- 
less spin-liquids occur in clean Mott insulators, where 
strong charge fluctuations and frustration prevent mag- 
netic ordering [7-9]. This experience makes it natural to 
ask whether or not one should expect a spin-liquid phase 
to form in (uncompensated) doped semiconductors near 
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FIG. 1. Two scenarios for the MIT in doped semiconduc- 
tors, (a) Conventional picture electron localization transition 
to Anderson-Mott insulator, which crosses over continuously 
to a pure Mott insulator (indicated by wiggly lines), (b) In 
the newly proposed scenario, the transition is to a spin-liquid 
with gapless fermionic spinon excitations, here the electri- 
cal MIT and spin/thermal MIT occur separately at n c i and 
n C 2 respectively. "+LM" indicates the fraction of randomly 
decoupled local spin-moments that inevitably accompany all 
phases. 

the MIT where charge fluctuations are strong, the system 
is at half-filling, and magnetic order is prevented by the 
random lattice structure, the competition between anti- 
ferromagnetic direct-exchange and random-sign RKKY 
exchange, and by quantum fluctuations. 

Nature of the Possible Spin-Liquid Phase(s) - The 
proposed spin-liquid phase is most conveniently de- 
scribed by formally dividing the electron into a bosonic 
U(l) rotor e l ° that carries the electron-charge, and a 
fermionic spinon f a that carries the electron-spin: Cj = 
e j6li /ia[13, 14]. This description allows extraneous un- 
physical states that must be removed by constraining 
n b.i ~ ^2 a flafia = 1 on each site, i. Here nu is the num- 
ber operator conjugate to 0j. The above decomposition 
has a U(l) gauge redundancy associated with 6i — > 



and/i 



1 fi, which manifests itself in the low-energy 



effective theory as an emergent U(l) gauge field, a(r,i) 
[15]. 

Decoupling the hopping term —Ujc\ s Cj. s = 
—tije l ( ei ~ 6 ^ fj 8 fj >8 in a mean field approximation 
and including gauge fluctuations gives the following 
effective action: S c g — J dr (Lb + Lf) 

Lb = T,ih + a ?) {uSij + v l3 y l (d T e 3 + 
-*«£«4- e<(fl< " 0<+Ow) 

Lf = E« fl [(dr a? M ) Sij - t id X%e-^\ f jt . (2) 

where x{j = (fi, s fj,s) and X % = (e i( - ei ~ 9 ^) are deter- 
mined self-consistently. Note that, due to the random 
placement of sites, in general ^2 s (fJ wm t> e spatially 
varying. Consequently, even at the mean-field level, the 
bosons will experience a random chemical potential; this 
changes the universality class of the Bose-Mott transition 
compared to the clean case (where (rib) = 1 f° r every site 
on both sides of the Mott transition). 

The metallic Fermi-liquid state corresponds to a 
superfluid-ordered phase of the bosonic rotors with 

( e *(0(x)+/ x a(y).dy)) ^ coex i stmg w j t h a diffusive 

Fermi-liquid of spinons. In this phase, the emergent 



gauge field is gapped by the condensate of charged rotors 
through the Anderson-Higgs mechanism, and the rotors 
and spinons are "glued" together into ordinary electrons. 

Eq. 2 also naturally describes a deconfined state in 
which the rotors form an insulator, while the spinons re- 
main diffusive. This results in an exotic charge-insulator 
with finite-density of states for gapless spin- 1/2 excita- 
tions. We suggest that this spin-liquid phase may occur 
near the MIT for doped semiconductors. In this sce- 
nario, shown in Fig. lb, the magnetic properties of the 
system change only gradually across the MIT at n c i , and 
are qualitatively identical in both the metal and insula- 
tor. In particular, we expect that one would still find 
a decoupled fraction of local-moments. As these local 
moments dominate the low-temperature thermodynam- 
ics and magnetic properties, the clearest signature of the 
spin-liquid is metallic thermal conductivity, k <~ T at low 
T[8]. While there has been extensive experimental analy- 
sis of the conductivity of doped semiconductors near the 
MIT, very little is known about thermal-transport. 

In the slave-rotor language, the formation of local 
moments comes from rare strong fluctuations in disor- 
der that locally bind the rotor and fermion back into 
a correlation-localized electron. We assume that the 
principal effect of correlated disorder among the rotor, 
spinon, and gauge-field sectors is to produce such local 
moments, and that the physics of the remaining non-local 
moment bulk can be well described by treating disorder 
separately in each sector. 

In the spin-liquid phase, the emergent gauge field is 
deconfined, and in clean systems its fluctuations lead to 
singular self-energies for the spinons resulting in non- 
Fermi-liquid behavior (2D) [16-23] or marginal Fermi- 
liquid behavior (3D) [24-26]. For the strongly disordered 
doped semiconductors, the inelastic scattering rate for 
the spinons from gauge fluctuations scales as r 9 _1 ~ T 
and is dominated by the clastic impurity scattering for 
low T[27]. Consequently, the low-energy properties of the 
disordered spinon Fermi-liquid will be largely unaltered 
by the emergent gauge field. Furthermore, the gauge field 
propagation is strongly damped by the diffusive spinons, 
leading to an u ~ q 2 scaling of gauge-excitations. This 
scaling implies that the gauge-field contribution to ther- 
modynamic quantities is sub-dominant compared to the 
spinon contribution. For example the gauge-field specific 
heat scales as C a ~ T 3 / 2 <c C sp inon ~ T. 

In 2D, a deconfined phase for the gauge-field requires 
the presence of extended, gapless fermionic excitations to 
suppress instanton configurations [28, 29]. In 3D, how- 
ever, a compact U(l) gauge-field may remain deconfined 
even without extended, gapless matter [28]. Therefore, 
in addition to the gapless, thermally-conducting spinon 
Fermi-surface state described above, it is also possible to 
form an insulating state where the charge degrees of free- 
dom are Mott localized and the spinons are Anderson- 
localized by disorder. Such a spinon Anderson insula- 
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tor is distinguished (in principle) from the conventional 
Anderson-Mott insulator by the presence of a gapless 
emergent U(l) gauge-field (though experimentally de- 
tecting the emergent gauge-field would be challenging). 

Generalized Phase Diagram - The MIT achieved by 
changing n, though experimentally relevant, is concep- 
tually complicated since disorder and interactions are 
simultaneously affected. It is conceptually simpler to 
consider a generalized phase-diagram where disorder 
strength W and interaction strength U can be separately 
adjusted, as in Fig. 2. Here we restrict our attention 
to 3D, half-filling, and non-nested Fermi-surfaces (which 
are not inherently unstable to magnetic ordering). Fur- 
thermore, we remain agnostic about the particular real- 
ization of disorder, with the expectation that such details 
will not alter the qualitative discussion that follows. 

We begin by considering various limiting cases. The 
(U = 0, W 7^ 0) limit is completely understood[12]: here 
a diffusive Fermi-liquid occurs up until a critical disorder 
strength beyond which all states near the Fermi-energy 
become localized leading to an Anderson insulator (AI). 
Each of these phases is known to be stable to infinitesi- 
mal interactions, and therefore extends at least to small 
U. The limit of (U ^ 0, W -> oo) is also straightfor- 
ward. Here the Anderson localized insulator at weak 
interactions crosses over continously to the Mott insu- 
lator of local-moments at strong-interactions. At T = 0, 
the local moments are magnetically ordered in either a 
random-singlet or spin-glass phase. 

Finally, the line (U ^ 0, W = 0) is also reasonably well 
understood [14], albeit with slightly less confidence. The 
clean Fermi-liquid survives up until some critical inter- 
action strength, beyond which it becomes a weak Mott- 
insulator with a spinon Fermi surface (SFS). For large U, 
the emergent gauge field undergoes a confinement transi- 
tion and anti-ferromagnetic order develops. Here again, 
each of the clean interacting phases is stable to infinitesi- 
mally small amounts of disorder and extends to finite W. 
The only distinction being that, for any (U ^ 0, W ^ 0), 
disorder creates a non-zero density of decoupled local mo- 
ments (indicated in Fig. 2 by "+LM"). 

These considerations greatly constrain the structure of 
the generalized phase-diagram. Each of the phases at 
the boundary are known to extend to finite values of W 
and U. Given the understanding of the boundaries of 
the phase diagram, the main issue here is not whether a 
strongly disordered fermionic spin-liquid could exist, but 
rather which particular path through the generic W and 
U phase-diagram is appropriate to tuning n in doped 
semiconductors. Fig. 2 depicts an extension of the well- 
understood outer boundary of the phase diagram to the 
interior. While one can conceive of many intermediate 
insulating phases at intermediate U and W, in the slave- 
rotor language, the only other natural candidate is the 
deconfincd spinon Anderson insulator described above. 

Thermal Conductivity - In the spin-liquid scenario, 
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FIG. 2. Schematic generalized phase-diagram as a function 
of disorder strength W and interaction strengths U measured 
with respect to the typical hopping t. 

the electrical MIT and thermal MIT occur separately: 
whereas electrical conductivity vanishes in the insulat- 
ing phase, thermal conductivity remains metallic, scal- 
ing as n sp ~ T at low temperatures. Since the ever- 
present concentration of local-moments dominates the 
low-temperature thermodynamical properties of the sys- 
tem (but contributes only weakly to transport), linear-T 
metallic conductivity is the clearest experimental signa- 
ture of the spin-liquid. 

While n sp <~ T at the lowest temperatures there will be 
Altshulcr-Aronov-type corrections to k from interactions 
and disorder: kaa ~ T 3 / 2 [30-32]. Also, one expects a 
large contribution, K p h from phonons: K p h ~ T 3 [39]. 
Therefore, to observe the metallic spinon-contribution, it 
may be necessary to work at very low temperature, and 
carefully subtract sub-dominant contributions. 

Quantum Critical (QC) Scaling - Despite extensive ex- 
perimental and theoretical effort, the quantum-critical 
(QC) behavior of electrical conductivity remains con- 
tentious and poorly understood. The existence of a spin- 
liquid phase would have important implications for how 
QC scaling should be extracted and interpreted. For 
T = and n > n c , the system is a Fermi-liquid obey- 
ing the Wiedemann- Franz law: k/LT = a (where the 
Lorenz number L is a constant). Since a vanishes at the 
transition while k/LT remains non-zero there must be 
a discontinuous jump in the T = electrical conductiv- 
ity at the MIT. In the slave-rotor description, this jump 
arises from the Ioffe and Larkin rule [33] that the electrical 
resistivity p equals the sum p = p b + pf of the resistivi- 
ties of the bosonic rotors pb and spinons pf respectively. 
Crossing the MIT at T = 0, the bosons transition from 
a superfluid with pb — to an insulator with p& = oo. In 
contrast, the fermionic contribution, pf evolves smoothly 
through the transition, implying a non-universal jump in 
the zero-temperature conductivity. Though superficially 
similar to Mott's early proposal[34], this jump in conduc- 
tivity is unrelated to the idea of a "minimum metallic 
conductivity" . 

Evidence against a discontinuous jump in conductiv- 
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FIG. 3. Quantum critical scaling (with z = 1)[40] of electrical 
conductivity a and linear-T coefficient of thermal conductiv- 
ity k/TL (L is the Lorenz number) near the MIT as a function 
of concentration and temperature. 

ity in Si:P comes mainly from pressure tuning studies[5] 
that show conductivity droping sharply but apparently 
continously to zero at the MIT. However, determining 
whether one is truly accessing the asymptotic behavior 
near the QC point is very challenging, and the proper 
interpretation of the conductivity scaling remains con- 
troversial and poorly understood [6]. For example, the 
pressure tuned experiments extract a conductivity scal- 
ing exponent v = 1/2 that is incompatible with general 
exponent inequalities for a metal-insulator transition [35], 
but could be explained as a thermally rounded version of 
the true T = conductivity jump. In the spin-liquid sce- 
nario, as we will argue below, issues with extrapolating 
to T — > and n —¥ n c are further exacerbated. 

Near the quantum critical point (QCP) (i.e. T « 0, 
and Sn = n — n c <C n c ), pb(T,S) obeys the quantum 
critical (QC) scaling for the disordered Bose-Hubbard 
model Mott transition as shown in Fig. 3a. In the high- 
temperature critical regime where T is the dominant per- 
turbation away from criticality, Pb(T) ~ T~ x / Z . At lower 
temperatures T < T* ~ (8ii) z l v , where Sn is the domi- 
nant perturbation from criticality pb(T) crosses over from 
the T-^ z to as superfluidity develops. Here z is the 
dynamical exponent for the disordered Bose-Mott tran- 
sition with Coulomb interactions [40]. 

The spinon contribution to the resistance tends to a 
constant at zero temperature, due to the elastic scat- 
tering from disorder. At finite temperature, there will 
also be non-constant contributions to pf. paa ~ VT, [36] 
and p p h ~ T~ 3-5 , from interactions and phonons respec- 
tively. The resulting electrical resistance, p = p& + pf, 
is depicted in Fig. 3b. for various Sn near the MIT. 
The main feature here, is the resistance upturn due to 
the nearly-critical fluctuating bosons, which quickly dis- 
appears below T < T* as p saturates to a non-universal 
constant set by pj. The corresponding T dependence of 



a = 1/p is shown in Fig. 3c. Notice the discontinous 
jump in the very-low temperature conductivity between 
Sn —> + and Sn = 0. As shown in Fig. 3d, this jump 
will be rounded at non-zero temperature, and could es- 
cape notice (consider, for example, if the lowest achiev- 
able temperature were indicated by the vertical dotted 
line in Fig. 3c). 

The spin-liquid scenario outlined here suggests a very 
different scheme for extracting the QC behavior of con- 
ductivity, than that for a conventional localization tran- 
sition. Here, one should include only data for which the 
resistance saturates to a nearly constant value set by the 
spinon contribution. In practice there is a minimum 
achievable value of temperature, T min . Consequently, 
this saturation region will disappear as the MIT is ap- 
proached when Sn < T^*. Beyond this point, extrapo- 
lations based on the curvature of a would fail to capture 
the true T — > behavior. 

The spin-liquid scenario will also complicate efforts to 
extract the critical scaling of a(n) near the MIT. This dif- 
ficulty is illustrated in Fig. 3d, which shows a(T = 0, n). 
As the concentration is decreased in the metal, the con- 
ductivity curves slowly towards an eventual localization 
transition at n C 2 (which may or may not occur). How- 
ever, in the present scenario, the Mott transition of the 
rotors intervenes at n c i > n c % before the spinons local- 
ize. In this case, extrapolations of QC scaling based on a 
conventional Anderson transition from the metallic side 
would be misleading. 

Discussion - In summary, we propose an alternative 
scenario to the Mott transition in doped semiconductors 
where the weakly insulating state is a spin-liquid with 
fermionic excitations. While such a transition has defi- 
nite consequences for the quantum critical scaling of con- 
ductivity near the MIT, such quantum critical behavior 
is notoriously difficult to determine. 

Other possible signatures of spin-liquid behavior in- 
clude sub-gap optical conductivity [37] in the insula- 
tor from gauge fluctuations and vanishing quasi-particle 
residue approaching the MIT (measurable by tunnel- 
ing on the metallic side). However, the former coexists 
with sub-gap conductivity from exciting weakly bound 
local moments, and the latter behavior will also be pro- 
duced by a soft Coulomb gap (which will develop at the 
MIT) [38]. Consequently, such probes are indirect, and 
would require a detailed quantitative comparison. 

Therefore, we suggest that the clearest test for spin- 
liquid behavior in doped semiconductors would come 
from a careful study of thermal transport near the MIT. 
A spinon Fermi-liquid would lead to k ~ T for low T, 
which, if observed, would strongly indicate the presence 
of gapless fermionic excitations. 
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Appendix — Disorder Averaged Properties of the 
Low— Energy Effective Slave— Rotor Theory 

In this appendix, we develop further, the properties 
of the low-energy effective field theory for the metal- 
insulator transition and spin-liquid phase within the 
slave-rotor theory. In principle, a full treatment of disor- 
der would need to account for correlations in the effects 
of local disorder fluctuations among the boson, fermion 
and the gauge field sectors. We assume that the princi- 
pal effect of such correlations is to produce local patches 
where the rotor and fermion are locally confined, i.e. to 
produce the dilute fraction of local-moments that are ob- 
served in each phase, but appear to decouple from the 
remainder of the system. We expect that the formation 
of local-moments accounts for the effects of the rare-long 
tails of the disorder distribution, and that the remaining 
connected component of the system is reasonably char- 
acterized by treating disorder separately in each sector. 

In the boson-sector, the superfluid-insulator transition 
will occur in the presence of a random potential, and the 
resulting insulating phase will be a mixture of glassy pud- 
dles of superfluid which do not percolate, coexisting with 
a Mott-localized bulk. In the spin-liquid scenario, the 
thermally conducting fluid of spinons (not including local 
moments) form a diffusive metal. As shown below, the 
phase transition in the boson sector is not affected by the 
presence of gapless fermions or gauge fluctuations, and 
is identical to that of the ordinary dirty Bose-Hubbard 
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model with random chemical potential and long-range 
Coulomb interactions [1]. 

A. Irrelevance of Spinon and Gauge Fluctuations on 
the Slave— Particle Boson— Mott Transition 



g' J s d+z d d qdws~ z G a s~ 2 \jb(u,q)\ 2 . Consequently g — > 
gs 2 ~ d flows to zero under RG and is irrelevant. Similarly, 
for z > 2, the term rescales as g J d d qdiuG a \jb(u!, q)| 2 — > 
g'J s d+z d d qdLos~ 2 G a s~' 2 \j b (uj,q)\ 2 , indicating 
g — »■ s 4 ~( d+z ^g, and the term is again irrelevant. 



Coupling to Fermions - Including gauge fluctuations 
gcncrically gives rise to a spinon-rotor density-density 
coupling of the form Ao#n/n&, where Srif is the devia- 
tion of rif from its average value. Here we give a simple 
scaling argument that such a coupling does not alter the 
critical behavior of the boson-sector near the Mott tran- 
sition. Consider integrating out the spinons. The leading 
order term in the effective action for the bosons will be 
of the form: A ^ iq (5n/(w, q)(5n/(-w, -q)) |n fc (cj, q)| 2 . 
In the transition to the spin-liquid, the spinon density- 
density correlator is diffusive and evolves smoothly across 
the transition: (5n,f(ui, q)5n/(— w, — q)) <~ \J\+D q 2 • 

After a momentum-shell renormalization-group (RG) 
step, integrating out modes with q e [A/s, A], with 
s > 1, one rescales q — > sq and u> — > s z uj to com- 
pare the new effective action to the original. Due to 
the random chemical potential provided by the spinon 
sector, the Boson Mott transition is in the same univer- 
sality class whether one tunes through the transition ei- 
ther by changing chemical potential or hopping strength 
(this would not be true without the presence of a random 
chemical potential, where the Bose-Mott transition takes 
place at fixed density per site, with the random potential 
however, the density per site is only fixed on average). 
For the chemical potential driven transition, under an 
RG step, the scaling part of the boson density rescales 
as ni,(r,t) s d+z ~ 1 / lJ nb{r,t) where v is the correlation 
length exponent. Equivalently, the fourier component 
rescales as nb(q,uj) — »■ s~ 1 /"rib(q, u>). 

For z < 2, the denominator of the diffusive fermion 
correlator is dominated by the \u>\ term, and the density 
coupling term scales as: Xg J d d qduq d ~ z \nb(q,Lj)\ 2 — > 
y s 2{d-i/v) j d d qdujq d - z \n b {q,w)\ 2 , indicating that the 
coupling constant A rescales as A — > s 2 ( 1 / l/ - £ \. Since 
I'd > 2, A is irrelevant. Similarly if z > 2, the diffusive 
fermion correlator scales like a constant under RG, and 
the coupling constant A — > s 2v ~ d ~ z \ is again irrelevant. 

Coupling to Gauge Fluctuations - The rotor-gauge 
field coupling generically takes the form J d d rdra tl j b , 
where % is the boson current. Integrating out the 
gauge-field at the RPA level generates a term of the 
form g J d d qdujG a \j b (uj,q)\ 2 , where G a = (|a(w, q)\ 2 ). 
Current-continuity requires that jb(u>,q) <~ ^rib(u),q), 
where n b is the total boson density which scales like 
n b ~ indicating that j&(w, q) — > s~ 1 j b (uj, q) in each 

RG step. As shown below, the gauge-field propagator 
scales like G a <~ (to + q 2 )^ 1 - 

For z < 2, under RG the gauge-fluctuation 
term rescales as g J d d qdu)G a \jb{u, q)| 2 — > 



B. RPA Effective Action for the Emergent Gauge 
Field in the Diffusive Spinon Metal 

Since the spinons and bosons have opposite charge un- 
der the emergent gauge field a, the gauge field couples 
to the currents as S a -j = J drd 3 r (j^ — jf) (here and 
throughout, we work in imaginary time). Therefore inte- 
grating out the spinon and boson fields within the RPA 
approximation, gives the following disorder-averaged ef- 
fective action for the gauge field: 

Seff PA) = ^a^u, q) K^(w,q) +Kf v (co,q) a v {u,q) 

u),q 

(3) 

where KfJ b are * ne disorder-averaged current-current 
correlators (equivalently linear-response kernels) for the 
fermions and bosons respectively. 

The temporal fluctuations of the gauge field arc 
screened by the compressible fermions and become mas- 
sive. At the critical point and in the insulating phase 
the Boson conductivity vanishes, and consequently the 
gauge field dynamics are determined by the fermion re- 
sponse. The disorder-averaged density-density part of 
the fermion electric response kernel is diffusive: 



2N(0)Dq 2 



" M + Dq 2 



(4) 



where, we have expanded the diffusion "constant" 
D(lj, q) (which generally has some w and q dependence, 
but does not vanish for q, uj — > outside of the disorder- 
localized phase) near w = = q, and dropped the 
irrelevant higher order terms. The other components 
of electric field response-function are related to K( 
by gauge invariancc and charge conservation: 



oo 



da " 



Kqq. Furthermore, 



,„ =pif J , and K\. 
we can identify 2N(0)D as the static uniform spinon- 
conductivity o~f. In addition, there is the usual diamag- 
netic response to fluctuating magnetic fields for 
Xciq 2 , where \d is the Landau diamagnetic susceptibil- 
ity of the spinons (whose average value is not altered by 
disorder) . 

Combining these considerations, and working in the 
Coulomb gauge (V • a = so that only the transverse 
gauge field a^ remains) gives the following RPA action 
for the gauge-field a: 



5 (RPA) = £ 



Xdq + 



(TfUJ 



u\ + Dq 2 



|a^( W , 9 )| 2 (5) 
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where D = is the diffusion constant (r is the disorder 
scattering time), and <jf is the spinon-conductivity. 



C. Inelastic Scattering of Spinons from Gauge 
Fluctuations 



Using the RPA expression for the gauge field propagator, one can find the leading (one-loop) self-energy for low- 
energy spinons near the Fermi-energy: 

S(iw,fc) = J duid 2 q±dq^D a (in,q) ^v F — ^j Gf(i(co — Cl),k + q) 



Gf{iu),\s) = - 



iu - £k + a^sgnw 
■ , . j til 



A,(fi,«) = 1 ff , oa (6) 



In the fermion Green's function £fc +9 <~ vp-gi + 0(<z 2 ), 
indicating that qii ~ f2. Consequently one may approxi- 
mate q 2 <C in the gauge field propagator, making D a a 
function of q± only. Furthermore, in this approximation, 

the current vertex (yp^-^ ~ v 2 F . 

Performing the q» integral then gives iV(0)sgn(u; — f2) 
independent of the disorder scattering time r. Consider- 
ing the case of lo > for definiteness, this limits the range 
of O integration from to u>. Since the dominant con- 
tributions come from fi <~ q\, one finds that the spinon 
self-energy due to diffusively-screened gauge fluctuations 



scales like: T,(iuj) <~ wlog(l/w) Continuing to real-time 
one finds the inelastic gauge field scattering rate scales 
like Im [T, r (uj)] ~ u, or equivalently T-~l lastic ~ T. At 
low-temperature, this inelastic scattering is clearly sub- 
dominant compared to the elastic impurity scattering. 
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